New constraint on the existence of the fi + — > e + 7 decay 
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The analysis of a combined dataset, totaling 3.6 x 10 14 stopped muons on target, in the search for 
the lepton flavour violating decay fi + —> e + 7 is presented. The data collected by the MEG experi- 
ment at the Paul Scherrer Institut show no excess of events compared to background expectations 
and yield a new upper limit on the branching ratio of this decay of 5.7 x 10~ 13 (90% confidence 
level). This represents a four times more stringent limit than the previous world best limit set by 
MEG. 

PACS numbers: 13.35.Bv; 11.30.Hv; 11.30.Pb; 12. 10. Dm 



The lepton flavour violating fi + — > c + 7 decay is pre- 
dicted to have an unobservable low rate within the Stan- 
dard Model of elementary particle physics (SM) , despite 
the existence of neutrino oscillations |l| . Conversely, the 
majority of new physics models Beyond SM (BSM), 
particularly in view of the recent measurements of a large 
6 13 at reactor and accelerator Q experiments, pre- 
dict measurable rates for this decay. An observation of 
the [i + — >• e + 7 decay would therefore represent an un- 
ambiguous sign of BSM physics, whereas improvements 
in the branching ratio upper limit constitute significant 
constraints on the parameter space, complementary to 
those obtainable at high energy colliders. 

The present best upper limit on the fi + —> e + 7 decay 
branching ratio B (B < 2.4 x 1CT 12 at 90% C.L.) was 
set by the MEG experiment [Io| with an analysis of the 
data taken in the years 2009-2010, for a total number of 
1.75 x 10 14 positive muons stopped on target. 

In this paper we present an updated analysis of the 
2009-2010 data sample, based on recently improved al- 
gorithms for the reconstruction of positrons and photons 



together with the analysis of the data sample collected in 
2011 with a beam intensity of 3 x 10 7 /i + /s, which corre- 
sponds to 1.85 x 10 14 stopped muons on target. Further- 
more the combined analysis of the full 2009-2011 statis- 
tics is presented. 

The signature of the signal event is given by a back- 
to-back, monoenergetic, time coincident photon-positron 
pair from the two body fi + —> c + 7 decay. In each 
event, positron and photon candidates are described by 
five observables: the photon and positron energies (i£ 7 , 
E c ), their relative directions (# e7 , <t>e-y) [HI and emission 
time (t cl ). Our analysis is based on a maximum likeli- 
hood technique applied in the analysis region defined by 
48 < E 1 < 58MeV, 50 < E c < 56MeV, |t e7 | < 0.7ns, 
\9&y\ < 50mrad and \4>e-y\ < 50mrad, which is described 
in detail in [lfj. We call "time sidebands" the regions 
in the variable space defined by 1< |i C7 | < 4 ns, U E 1 - 
sideband" that defined by 40 < E 1 < 48 MeV and "an- 
gle sidebands" those defined by 50 < |0 C7 | < 150mrad or 
50 < |0 C7 | < 150mrad. 

The background has two components, one coming 
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from the Radiative Muon Decay fi + — > e + vD-f (RMD) 
and one from the accidental superposition of energetic 
positrons from the standard muon Michel decay with 
photons from RMD, positron-electron annihilation-in- 
flight or bremsstrahlung. At the MEG data taking rate, 
93% of events with E 1 > 48 MeV are from the Acciden- 
tal background (ACC). 

The MEG detector is described in detail elsewhere 12 1 . 
It is comprised of a positron spectrometer formed by a 
set of Drift CHambers (DCH) and scintillation Timing 
Counters (TC), located inside a superconducting solenoid 
with a gradient magnetic field along the beam axis, and 
a photon detector, located outside of the solenoid, made 
up of a homogeneous volume (900 £) of Liquid Xenon 
(LXe) viewed by 846 UV-sensitive photomultiplier tubes 
(PMTs) submerged in the liquid. 

The MEG detector response, resolutions and stability 
are constantly monitored and calibrated by means of a 
multi-element calibration system 1^, 12j, which was re- 
cently complemented with a new method based on the 
Mott scattering of a monochromatic positron beam. Two 
important hardware improvements were introduced in 
2011, one in the 7r~p — > 7r°n Charge EXchange reaction 
(CEX) calibration by replacing the Nal detector, used to 
define the back-to-back coincidence with the LXe detec- 
tor when observing the two photons from ir° decay, with 
a higher resolution BGO array detector, the other in the 
optical survey technique for the DCH by using a laser 
tracker and prismatic corner cube reflectors mounted on 
the DCH modules. 

In 2011 the DAQ efficiency, being the product of DAQ 
live time fraction and trigger efficiency to select signal 
event, was at the level of 96%. An improvement from 
72% in 2010 is due to a new multiple buffer read out 
scheme, which guarantees a 99% live time fraction with 
a less constrained event selection, thus allowing a 97% 
signal trigger efficiency [l3[ • 

The positron track is reconstructed by combining its 
measured positions at each DCH layer (hit) in the spec- 
trometer. Longitudinal (along the muon beam direction) 
z-positions are derived from signals induced on the seg- 
mented DCH cathodes. Electromagnetic noise has been 
the main source of degradation of the z-resolution dur- 
ing all the data-taking periods and a new reconstruction 
algorithm based on a fast Fourier transform filtering tech- 
nique has been applied to mitigate such effect, yielding an 
up to ~10% improvement in angular resolution. Internal 
alignment of the DCH layers is obtained by tracking cos- 
mic ray muons with the magnet off and using an optical 
survey, as described in }10| . 

The positron kinematic variables are extracte d by 
means of a Kalman filter track fitting technique 1J, [l_5| . 
This algorithm has been completely revised for this anal- 
ysis to include a better model for the hits and the track it- 
self, based on the GEANE package 1|| 17 1. An improved 
model for the detector material, accounting for multiple 



scattering and energy loss as well as a detailed map of the 
non-uniform magnetic field, measured with a 0.2% pre- 
cision has also been included. The fitted positron track 
is propagated to the TC allowing an iterative refinement 
of the hits with the positron time measurement. 

The track fit yields a covariance matrix for the parame- 
ters, resulting in very good agreement with the measured 
resolutions, which are extracted from a sample of tracks 
with two full turns in the DCH, by comparing the track 
parameters determined independently for each turn at 
an intermediate plane. Consequently a per-track error is 
determined which allows us to follow the variable DCH 
performance during the data-taking period and is taken 
into account in the maximum likelihood analysis. 

The average hit multiplicity for a track is about 10 
and only tracks with at least 7 hits and either one or two 
turns in the spectrometer are retained for the analysis. 
Additional quality requirements on the \ 2 fit value and 
parameter uncertainties are applied to select only one 
positron per event. The overall improvement in positron 
reconstruction with respect to the previous data analysis 
is clearly visible in FigQ] (top) where the reconstructed 
positron energy near the kinematic edge of the Michel 
decay spectrum shows a reduced tail. The energy reso- 
lution is evaluated by fitting the kinematic edge and it is 
well described by the sum of three Gaussian curves with 
a resolution of <je — 305 keV for the core component 
(85%). 

The resolution in the azimuthal angle of the positron 
when it leaves the target, <j> c , has a dependence on <j} e 
itself with a minimum at (f> e = 0, where it is measured 
by the two-turn method after correcting for all known 
correlations to be — 7.5 (7.0) mrad [l8(. Similarly, the 
resolution in the polar angle, 9 e , is measured to be ag B 
— 10.6 (10.0) mrad. The decay vertex coordinates and 
the positron direction at the vertex are determined by 
extrapolating the reconstructed track back to the target. 
The resolutions on the decay vertex coordinates are also 
determined by the two-turn method and are described 
by a Gaussian curve with a z = 1.9 (1.5) mm and, in the 
vertical direction, by the sum of two Gaussian curves 
with u y — 1.3 (1.2) mm for the core component (85%). 

The LXe detector uses the xenon scintillation light to 
measure the total energy released by the photon, the 
position and time of its first interaction. The three- 
dimensional photon conversion point is reconstructed by 
using the distribution of the number of scintillation pho- 
tons detected by the PMTs near the incident position. 
The photon direction is defined by the line connecting 
the decay vertex to the photon conversion point in the 
LXe detector. The reconstruction of the photon energy 
is based on the sum of the scintillation light detected 
by all PMTs. Monochromatic 55 MeV photons from 7r° 
decays are used to determine the absolute energy scale. 
The photon conversion time is reconstructed by combin- 
ing the leading edge times of the PMT waveforms. 
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It is important to identify and unfold photon pile-up 
events at high muon rates since at 3 x 10 7 /i + /s beam 
rate, around 15% of triggered events suffer from pile- 
up. For the previously published analyses, photon pile- 
up events were identified topologically by the pattern of 
PMT light distribution and temporally by the leading 
edge time distribution in the time reconstruction, with- 
out the use of detailed waveform information. In addi- 
tion to these methods, a new algorithm, analyzing wave- 
forms after summing up all channels at the end of the 
full chain of photon reconstruction, was developed. It 
enables the efficient identification and removal of pile- 
up photons by using template waveforms. Consequently 
the charge integration window for the energy estimate is 
re- adjusted, resulting in a better energy reconstruction. 
This improvement is shown in Fig. [T] (bottom) where we 
compare the photon energy spectra obtained with differ- 
ent pile- up elimination algorithms. The reduced tail in 
the high energy region is clearly visible. The efficiency of 
photon reconstruction is improved from 59% to 63% due 
to the new algorithm. 

The performance of the position reconstruction is eval- 
uated by a Monte Carlo simulation, resulting in resolu- 
tions of 5 mm on the photon entrance face, and 6 mm 
along the radial depth from the entrance face. This is 
validated in CEX runs by placing lead slit collimators in 
front of the LXe detector. The timing and the energy 
resolutions are evaluated using two simultaneous back- 
to-back photons from w° decay. The LXe timing reso- 
lution is found to be at = 67 ps at the signal energy. 
The position-dependent energy resolutions are measured 
in the CEX data and the average energy resolution ex- 
tracted from a Gaussian fit to the high energy side of 
the spectrum is evaluated to be 1.7% (1.9%) and 2.4% 
(2.4%) for radial depths larger and smaller than 2 cm re- 
spectively. These position-dependent energy resolutions 
are incorporated into the likelihood analysis. 

The resolutions of the relative directions (6* C7 , 4> c -y) 
are derived by combining the relevant resolutions of 
positrons and photons discussed above. The results are 
16.2 (15.7) mrad for 6 C1 and 8.9 (9.0) mrad for <fi ei after 
correcting for all known correlations. The relative time 
t CJ is derived from the time measurements, one in the 
LXe detector and the other in the TC, after correcting 
for the lengths of the particle flight-paths. The associated 
resolutions at the signal energy are 127 (135) ps, evalu- 
ated from the RMD peak observed in the i? 7 -sideband; a 
small correction takes into account the ^-dependence 
measured in the CEX calibration runs. The position 
of the RMD-peak corresponding to t cl = was moni- 
tored constantly during the physics data-taking period 
and found to be stable to within 15 ps. 

A blind analysis procedure is applied only to the new 
dataset in 2011 by masking a region of 48 < E 1 < 58 MeV 
and \t cl \ < 1 ns until the Probability Density Functions 
(PDFs) for the likelihood function are finalized. For all 
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FIG. 1: (Top) The Michel positron spectrum in the same 
dataset in the time sidebands with the old (black dashed 
line) and the new (red solid line) track reconstruction code. 
(Bottom) The photon background spectra from the time side- 
bands of the muon data in 2010 with different pile-up elim- 
ination algorithms. Black dots: no pile-up elimination; blue 
dot-dashed: previous algorithm; red solid: new algorithm. 



the datasets including 2009-2010 data the background 
studies and the extraction of the PDFs are carried out 
in the time and angle sidebands. The maximum likeli- 
hood fit is performed in order to estimate the number 
of signal, RMD and ACC events in the analysis region. 
The definition of the likelihood function is described in 
detail in [l(| • All PDFs as a function of the observables 
are extracted from the data. Different resolutions and 
correlations are used in the PDFs on an event- by- event 
basis. The dependence on the photon interaction posi- 
tion and the quality of the positron tracking has already 
been incorporated into the previous analysis, while in the 
new analysis a per-event error matrix for the positron ob- 
servables, estimated by the new Kalman filter, has been 
introduced into the PDFs. The sensitivity is improved 
by about 10% in the new analysis with the positron per- 
event error matrix. An analysis with constant PDFs is 
also performed as a crosscheck, showing consistent re- 
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suits. The confidence interval for the number of signal 
events is calculated by a frequentist method with a profile 
likelihood-ratio ordering ID, 2(| , where the numbers 
of RMD and ACC events are treated as nuisance param- 
eters. 

To translate the estimated number of signal events into 
a signal branching ratio two independent normalization 
methods are used, either counting the number of Michel 
positrons selected with a dedicated trigger or the num- 
ber of RMD events observed in the muon data. Their 
combination leads to a 4% uncertainty in the branching 
ratio estimate. The increased reconstruction efficiency of 
the new algorithms results in a 14% larger data sample 
for the n + — > e + 7 search, as estimated with both nor- 
malization methods, with both positron and photon new 
algorithms contributing equally. 

The systematic uncertainties on the PDF parameters 
and on the normalization are taken into account in the 
calculation of the confidence intervals by fluctuating the 
PDFs by the amount of the uncertainties. In total they 
produce a 1% effect on the observed upper limit, with 
the majority of the contribution coming from the angular 
PDFs. 

The sensitivity (1S90) is estimated as the median of the 
distribution of the branching ratio upper limits at 90% 
C.L., calculated over an ensemble of pseudo-experiments, 
randomly generated according to the PDFs based on a 
null signal hypothesis, with the rates of ACC and RMD 
evaluated from the sidebands. The sensitivities have been 
so evaluated for the 2009-2010 combined data, the 2011 
data alone and the 2009-2011 combined data sample, and 
are reported in Table [U Likelihood analyses are also per- 
formed in fictitious analysis regions in both the time- and 
angle-sidebands, getting upper limits all in good agree- 
ment with the iStjo's. 

Figure [2] shows the event distributions in the (E^^E^)- 
and (cos C7 , t C 7)-planes for the combined 2009-2011 
dataset, where Q ei is the opening angle between positron 
and photon, together with the contours of the averaged 
signal PDFs. 

The observed profile likelihood ratios as a function of 
the branching ratio are shown in Fig. [31 The best B es- 
timates, upper limits at 90% C.L. (#90) and S90 for the 
combined 2009-2010 dataset, the 2011 data alone and the 
total 2009-2011 dataset are listed in TableQ] The B 90 for 
the latter is 5.7 x 10~ 13 . As a quality check the maximum 
likelihood fit is repeated for the 2009-2011 dataset omit- 
ting the constraint on the number of background events. 
We obtain A/rmd = 163 ± 32 and iV A cc = 2411 ± 57, 
in good agreement with the expectations estimated from 
Ej and time sidebands, (Nrmd) = 169 ± 17 and (-/Vacc) 
= 2415 ± 25. 

The reanalysis of the 2009-2010 dataset with new al- 
gorithms has led to variations in the values of the ob- 
servables which are much smaller than the detector res- 
olutions. The events observed with the highest signal- 
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FIG. 2: Event distributions for the combined 2009-2011 
dataset in the (E B ,E^)- and (cos Q el , t 07 )-planes. In the 
top (bottom) panel, a selection of \tej\ < 0.244ns and 
cos C7 < —0.9996 with 90% efficiency for each variable 
(52.4 < E c < 55 MeV and 51 < E 1 < 55.5 MeV with 90% 
and 74% efficiencies for E e and E ly respectively) is applied. 
The signal PDF contours (1, 1.64 and 2 a) are also shown. 



likelihood in the previous analysis of the 2009-2010 
dataset have also moved in the new analysis within the 
expected fluctuation and mostly have had their signal- 
likelihood slightly reduced. These small variations induce 
a change in B&t and Bqq for the same dataset. We have 
compared $90 's obtained with the new and old analyses 
for the same sample of simulated experiments and found 
that a change of Bgo equal to or larger than what we 
observe in the 2009-2010 dataset has a 31% probability 
of occurring. The upper limit obtained from the 2011 
data only is more stringent than Sgo- This is, however, 
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TABLE I: Best fit values (Bat's), branching ratios (Bgo) and 
sensitivities (<Sgo) 



program is underway aiming at a sensitivity improvement 
of a further order of magnitude 21 1. 



Dataset 


S flt x 10 12 


B 90 x 10 12 


Sgo x 10 12 


2009-2010 


0.09 


1.3 


1.3 


2011 


-0.35 


0.67 


1.1 


2009-2011 


-0.06 


0.57 


0.77 
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FIG. 3: Observed profile likelihood ratios (A p ) as a function 
of the branching ratio for the 2009-2010 combined data, the 
2011 data alone and the combined 2009-2011 data sample. 



not considered unusual, since the probability to have Bqq 
equal or smaller than that observed in the 2011 data is 
calculated to be 24% with a sample of simulated experi- 
ments. 

In conclusion the MEG experiment has so established 
the most stringent upper limit to date on the branching 
ratio of the /i+ e+7 decay, B < 5.7 x 10~ 13 at 90% 
C.L. using data collected between 2009 and 2011, which 
improves the previous best upper limit by a factor of four. 
Further data have also been acquired in 2012 with an 
additional three-month run scheduled for 2013; the final 
number of stopped muons is expected to be almost twice 
that of the sample analyzed so far. Currently an upgrade 
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